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Abstract. The model of holographic dark energy (HDE) with massive neutrinos and/or dark 
radiation is investigated in detail. The background and perturbation evolutions in the HDE 
model are calculated. We employ the PPE approach to overcome the gravity instability 
difficulty (perturbation divergence of dark energy) led by the equation-of-state parameter w 
evolving across the phantom divide re = —1 in the HDE model with c < 1. We thus derive the 
evolutions of density perturbations of various components and metric fluctuations in the HDE 
model. The impacts of massive neutrino and dark radiation on the CMB anisotropy power 
spectrum and the matter power spectrum in the HDE scenario are discussed. Eurthermore, 
we constrain the models of HDE with massive neutrinos and/or dark radiation by using the 
latest measurements of expansion history and growth of structure, including the Planck CMB 
temperature data, the baryon acoustic oscillation data, the JLA supernova data, the Hubble 
constant direct measurement, the cosmic shear data of weak leasing, the Planck CMB leasing 
data, and the redshift space distortions data. We find that < 0.186 eV (95% CL) and 

IVeff = 3.75;0;1 in the HDE model from the constraints of these data. 
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1 Introduction 

Detection of neutrino oscillations has indicated that neutrinos have masses. However, the 
neutrino oscillation experiments cannot measure the absolute masses of neutrinos, but can 
only measure the squared mass differences between the neutrino mass eigenstates. For exam¬ 
ple, the solar and reactor experiments observed ~ 8 x 10“^ eV^, and the atmospheric 

and accelerator beam experiments observed Am ‘^2 — 3 x 10“^ eV^. As a complementary 
to laboratory experiments, current cosmological data have been used to get information on 
the absolute scale of neutrino masses. Actually, the current available cosmological data have 
been providing tight limits on the total mass of neutrinos. See Refs. [1-3] for reviews of 
neutrino cosmology. 

The cosmic microwave background (CMB) observations have been used to constrain the 
neutrino mass and the extra relativistic degrees of freedom (sometimes referred to as “dark 
radiation”) [4-6]. In the base 6-parameter ACDM model, a normal mass hierarchy with 
Y2 ~ 0.06 eV (dominated by the heaviest neutrino mass eigenstate) is assumed. When 
the base ACDM model is extended to allow for larger neutrino masses, one often assumes 
three species of degenerate massive neutrinos, neglecting the small differences between mass 
eigenstates. The Planck CMB temperature power spectrum provided the tight limits on 
the total mass of active neutrinos, and the effective number of relativistic species, 

Aeff [7]. For example, the Planck-|-WP-|-highL data combination (here, WP denotes the 
WMAP 9-year polarization data, and highL denotes the ACT and SPT temperature data) 
gives the 95% confidence level (CL) limits: < 0.66 eV for the case of no extra relics 

(Aeff = 3.046) and Aeff = 3.36;°;®® for the case of minimal-mass normal hierarchy for the 
neutrino masses (only one massive eigenstate with rriu = 0.06 eV). Note here that, unless 
otherwise specihed, Planck data in this paper refer to the Planck 2013 release. Late-time 
geometric measurements can be used to help reduce some geometric degeneracies and thus 
improve constraints. Therefore, the baryon acoustic oscillation (BAO) data are very useful 
in the parameter estimation. The Planck-|-WP-|-highL-|-BAO data combination changes the 
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above limits to: Y1 < 0-23 eV for the case of no extra relics and A^eff = 3.301 q 5 ^ for the 
case of minimal-mass normal hierarchy model. ^ 

Since the Planck data are in tension with the weak leasing measurements and the abun¬ 
dance of rich clusters, recently there has been significant interest in larger neutrino masses 
because neutrino free streaming provides a possible way to lower the late-time fluctuation 
amplitude cg and thereby could reconcile the above tensions. For the relevant discussions, 
see, e.g.. Refs. [9-17]. But it was also argued by some authors [18] that the larger neutrino 
masses offer only a marginal improvement compared to the base ACDM model. Anyway, 
the use of the large-scale structure measurements could help constrain the neutrino mass 
significantly. 

The cases of neutrinos and dark radiation in a dynamical dark energy model have been 
discussed [19]. In Ref. [19], the ACDM-l-neutrino/dark radiation models are extended by 
replacing the cosmological constant with the dynamical dark energy with constant w. The 
corresponding u)CDM-based models are also constrained by using the recent observational 
data in Ref. [19]. In this paper, we will further extend the relevant discussion. We will 
consider the cases of neutrinos and/or dark radiation in the holographic dark energy model. 

The holographic dark energy (HDE) model [20] originates from the holographic principle 
of quantum gravity, and so it remains significant interest in cosmology. It is expected that 
theoretical and phenomenological studies on HDE might provide important clues for the 
bottom-up exploration of a full quantum theory of gravity. There have been numerous studies 
on the theoretical implications and observational constraints of the HDE model [21-42]. In 
particular, the neutrino mass in the HDE model has been constrained in the light of the 
WMAP 7-year data in Ref. [39]. In this paper, we will make a more sophisticated analysis; 
we will discuss the constraints on the neutrino mass and/or the extra relativistic degrees 
of freedom in the HDE model by using the Planck CMB data plus other measurements of 
expansion history and growth of structure. 

We will consider three models in this paper, i.e., (i) the HDE-|-X]^i^ model, (ii) the 
HDE-|-Aefr model, and (hi) the HDE-|-^ mi,-|-Aefr model. We will constrain the models by 
using the latest observational data. The basic data combination we use is the CMB-I-BAO. 
Purthermore, we combine other geometric measurements, i.e., the type la supernovae (SN) 
JLA data and the direct determination of the Hubble constant Hq. Einally, we also consider 
to further combine the measurements of growth of structure, i.e., the weak lensing (WL) data 
and the redshift space distortions (RSD) data. 

The paper is organized as follows. In Sec. 2, we briefly describe the model of holographic 
dark energy with massive neutrinos and dark radiation. Por the cases of c < 1 (note that here 
c is not the speed of light, but the parameter of HDE) in the HDE model, the dark-energy 
equation-of-state parameter (EOS) w crosses the phantom divide w = —1. Usually, the 
perturbation instability appears when w crosses —1. In order to overcome this difficulty, we 
employ the “parametrized post-Eriedmann” (PPP) approach [43] to treat the perturbations 
in dark energy. We will illustrate the density perturbations of various components and the 
metric fluctuations in the HDE model. Also, we will discuss the impacts of massive neutrinos 
and dark radiation on the CMB anisotropy spectrum and the matter power spectrum. In 

^During the completion of this paper, the Planck 2015 results appeared on arXiv [8] (but the new data 
have not been released currently). According to the latest Planck 2015 data, it is found that 'Y^rrii, < 0.72 
eV from Planck TT-|-lowP and < 0.21 eV from Planck TT+lowP-|-BAO; = 3.13 ± 0.32 from 

Planck TT+lowP and A7ff = 3.15±0.23 from Planck TT-|-lowP-|-BAO. Here “lowP” denotes the Planck low-£ 
temperature-polarization data. 
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Sec. 3, we use the latest observational data to constrain the models. We will focus on the 
constraint results of c, flm, Hq, as, and iVeff- Conclusion is given in Sec. 4. 

2 A brief description for the model of holographic dark energy with mas¬ 
sive neutrinos and dark radiation 

The cosmological constant suffers from the severe theoretical challenge, i.e., it cannot be 
understood why the theoretical value of A is greater than the observational value by many 
orders of magnitude. The cosmological constant is equivalent to the vacuum energy density, 
and thus its value is determined by the sum of the zero-point energy of each mode of all the 
quantum fields. Thus we have the vacuum energy density pA — , where /cmax 

is the imposed momentum ultraviolet (UV) cutoff. When taking the UV cutoff to be the 
Planck scale (10^® GeV), where the quantum field theory in a classical spacetime metric is 
expected to breakdown, the vacuum energy density would exceed the critical density by some 
120 orders of magnitude. 

Obviously, one should not calculate the value of A in the context without gravity. In 
the absence of a full theory of quantum gravity, one has to explore the potential solution of 
the cosmological constant problem through combining the holographic principle of quantum 
gravity with the effective theory of quantum fields. That is the origination of the HDE model. 

Considering gravity in a quantum field system, the conventional local quantum field 
theory would break down due to the too many degrees of freedom that would cause the 
formation of a black hole. However, once the holographic principle is considered, the number 
of degrees of freedom can be reduced. 

One could put an energy bound on the vacuum energy density, paLP' < Mp^L, where 
Mpi is the reduced Planck mass, which implies that the total energy in a spatial region with 
size L should not exceed the mass of a black hole with the same size [44] . The largest length 
size compatible with this bound is the infrared (IR) cutoff size of this effective quantum field 
theory. Evidently, the holographic principle gives rise to a dark energy model based on the 
effective quantum field theory with a UV/IR duality. Erom the UV/IR correspondence, the 
UV problem of dark energy can be converted into an IR problem. A given IR scale can 
saturate that bound, and thus one can write the dark energy density as [20] 

PA = Sc^MpjL (2.1) 

where c is a dimensionless parameter characterizing all of the uncertainties of the theory. 
This indicates that the UV cutoff of the theory would not be fixed but run with the evolution 
of the IR cutoff, i.e., femax oc 

The HDE model chooses the event horizon of the universe 

RMt) = ait) / — (2.2) 

Jt a{t') 

as the IR cutoff L of the theory, explaining the fine-tuning problem and the coincidence 
problem at the same time in some degree [20] . Actually, it is clear to see that the holographic 
dark energy is essentially a holographic vacuum energy. However, this holographic vacuum 
energy does not behave like a usual vacuum energy, owing to the fact that its equation of 
state parameter w is not equal to —1. 

The HDE is a dynamical dark energy. The parameter c plays a crucial role in deter¬ 
mining the evolution of the HDE. At the early times (t —>• 0 or z —)• oo), the EOS of HDE 
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w —)> —1/3. The value of c determines the subsequent evolution of w. In the infinite future 
(t —>■ oo or 2 ; —)■ —1), we have w —>• —1/3 — 2/(3c). Thus, when c < 1, we find that w of HDE 
crosses —1 during the cosmological evolution. 


2.1 Background and perturbation evolutions of holographic dark energy 

Consider a spatially flat universe with HDE, matter, and radiation. The Friedmann equation 
reads 

3Mpji7^ = PA Pc + Pb + P'y Pu, (2-3) 

where pA, pc, Pb, P^y, and pi, stand for the energy densities of HDE, cold dark matter (CDM), 
baryons, photons, and neutrinos. The CDM and baryons are non-relativistic matter, and 
the photons are of course radiation. The neutrinos contribute to the radiation density at the 
early times, but they behave as matter after the non-relativistic transition. Actually, there 
may be some extra radiation, usually referred to as “dark radiation”. Their contribution to 
the total radiation content can be parametrized in terms of the effective number of neutrinos, 
Aeff- So the energy density of total radiation can be given by the following relation 


Pr — 


7 

CO 


TT 

8 

viiy J 


P'y: 


(2.4) 


where we see that pr has been normalized by the energy density of photons due to the fact 
that its present-day value is measured from the CMB temperature. 

Since the HDE model has been investigated extensively, here we directly give the back¬ 
ground evolution equations that are coupled differential equations governing the functions 
Da( 2 ;) and E{z). The equations are [39, 40]: 


1 dE{z) _ VIa{z) ( 1 / o ^ ^ , ^a(^) - ^r(^) 

E{z) dz l-|-z\c 2VLa{z) 

d^A{z) 2 Da( 2:)(1 - Da( 2 )) /I . 1 - ^Aiz) + nriz)\ 

=- 2 ( 1 -OaW) )• 

where E{z) = H{z)/Hq is the dimensionless Hubble expansion rate, and ^,.( 2 ;) 
z)'^/E{z)‘^. The initial condition for solving the above differential equations are: 
and Da( 0 ) = 1 — — D^. 

The EOS of HDE is given by the following relation [20] 

^(^) = ^ - (2-7) 

Thus, once the function Vt.A{z) is derived as a solution to the differential equations (2.5) and 
( 2 . 6 ), the evolution of w{z) is also entirely determined. 

The linear metric and matter density perturbations of the HDE model can be calculated 
by using the formalism of Ma and Bertschinger [45]. But when treating the dark energy 
perturbations, one should be very careful about the divergence problem at the w = —1 
crossing. We follow the WMAP and Planck collaborations to deal with this issue by using 
the PPF approach [43]. 

The perturbations of dark energy can be described by four variables: density fluctuation 
(5pdej velocity Ude, pressure fluctuation (5pde) and anisotropic stress Hde- For these quantities, 
the evolution equations of (5/9de and Ude are determined by continuity and Navier-Stokes 


( 2 . 6 ) 
= Dr(l -h 

: E(0) = 1 
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Figure 1. The evolutions of matter and metric perturbations in the HDE model at k= 0.01 Mpc“^, 
0.1Mpc“^, and 1.0Mpc“^. Here the matter density perturbations are calculated in the synchronous 
gauge, and the metric perturbations $ and dt are the gauge-invariant variables. We fix ^ rrii, = 0.06 
eV and A^eff = 3.046, and other parameters are fixed to be the best-fit values from Planck. 


equations, and 11^6 vanishes for the linear perturbations. In usual fluid approach, to complete 
the system, one need to specify the relationship between 6pde and 6pde by defining the rest- 
frame sound speed for dark energy. But the PPF method abandons this condition. The 
PPF approach replaces the condition of pressure perturbation with a direct relationship 
between the momentum densities of dark energy and other components on large scales, which 
determines the velocity Ude- Actually, once Vde is determined, Spde follows by the momentum 
conservation, and in this method no divergence will appear when w crosses the phantom 
divide. The details about the PPF description in the HDE model can be found in Sec. 2.2 
of Ref. [39]. For more details for the PPF code, we refer the reader to Ref. [43]. Note also 
that the PPF framework for interacting dark energy has been established recently [46, 47]. 

Figure 1 shows a concrete example of the numerical calculation for the evolution of mat¬ 
ter density perturbations and metric fluctuations in the HDE model with c = 0.8. To make 
a comparison for different scales, we choose three scales, i.e., k = 0.01, 0.1, and 1.0 Mpc“^, 
to illustrate. In this example, other parameters are fixed as the values in the 6-parameter 
ACDM model fitting to the Planck data. In particular, we note that the normal hierarchy 
minimal mass assumption is made for the neutrino masses, i.e., '^rriu = 0.06 eV, and Agfr 
is fixed to be 3.046. We show the evolution of the density perturbations of HDE, CDM, 
baryons, photons, and massless neutrinos in the synchronous gauge. Outside the horizon, 
the behavior of the density perturbations is strongly gauge-dependent. In the synchronous 
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Figure 2. The CMB anisotropy power spectrum Cf'^ and the matter power spectrum P{k) in the 
HDE model with c = 0.8. In the upper panels, the parameter '^rrii, is varied and other parameters 
are fixed; we choose '^rriv = 0, 0.5 eV, and 1.0 eV. In the lower panels, the parameter is varied 
and other parameters are fixed; we choose AlVeg = 0, 0.5, and 1.0. 


gauge, before the horizon crossing, all the density contrasts grow; after the horizon crossing, 
they come into causal contact and become nearly independent of the coordinate choices. For 
the HDE, we find that after the horizon crossing, the density perturbations in HDE become 
constant on large scales (e.g., k = 0.01 Mpc“^) and decrease quickly on small scales (e.g., 
A; = 0.1 and 1.0 Mpc“^). Hence, the HDE would not cluster significantly on the sub-horizon 
size. For all the cases, density perturbation of HDE is always smaller than that of CDM 
by several orders of magnitude, and is nearly in the same order of magnitude with that of 
radiation, and so the HDE perturbation almost does not affect the evolution of the CDM 
perturbation. The evolutions of the metric fluctuations and T are also shown in this figure. 

2.2 Impacts of massive neutrinos and dark radiation on CMB anisotropy spec¬ 
trum and matter power spectrum 

Neutrinos do not interact much at all with baryons for z <C 10^*^ (after thermal decoupling). 
As a result, neutrinos are treated as free-streaming particles. Since dark radiation behaves like 
massless neutrinos, it is also treated as free-streaming. They can affect the CMB anisotropy 
spectrum and matter power spectrum, thus providing a potential way to constrain them 
through CMB and large-scale structure (LSS) observations. 
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Massive neutrinos (with the total mass in the range from 10“^ eV to 1 eV) should be 
viewed as radiation at the time of equality and as non-relativistic matter today. The change 
in the redshift of equality Zeq affects the position and amplitude of the peaks. The change in 
the non-relativistic matter density at late times leads to the two effects: impact the angular 
diameter distance to the last scattering surface, Da{z^), that controls the overall position 
of CMB anisotropy spectrum features, and affect the slope of the \ow-i tail of the CMB 
spectrum through the late Integrated Sachs-Wolfe (ISW) effect. 

Neutrinos with masses can affect the matter power spectrum. On large scales [k < fcnr, 
with fenr — h Mpc“^), the neutrino free-streaming can be ignored, and 

so the neutrino perturbations are nearly indistinguishable from the CDM ones. Thus, the 
neutrino mass almost does not impact on the large-scale power spectrum. But the small- 
scale matter power spectrum P{k) with k > knr is reduced by the massive neutrinos. Due to 
the free-streaming of neutrinos, massive neutrinos do not cluster on small scales. Also, the 
growth rate of CDM perturbations is reduced via the absence of gravitational back-reaction 
effects of free-streaming neutrinos. 

The dark radiation (parametrized by AAgfr = Agfr ~ 3.046) affects the CMB power 
spectrum in several ways. First, the extra radiation density increases the early expansion 
rate, so the angular size of the acoustic scale, 0 * = Ts/D^, is reduced, which determines the 
peak positions. Second, the extra dark radiation density delays the time of equality, and so 
enhances the first and second peaks due to the early ISW effect. Third, as a free-streaming 
fluid, the dark radiation has a non-negligible anisotropic stress, and this changes the metric 
fluctuations during the radiation era and thus the temperature fluctuations on scales i > 130, 
because these scales enter the horizon during the radiation-dominated era. Finally, the extra 
radiation density increases the expansion rate and thus increases the diffusion length relative 
to the sound horizon, which enhances the Silk damping of the small-scale anisotropy. The 
dark radiation affects the matter power spectrum is also through the free-streaming, thus 
only the small-scale powers are suppressed by larger Agff- 

Figure 2 shows an example of how massive neutrinos and dark radiation impact on the 
CMB anisotropy power spectrum and the matter power spectrum P{k) in the HDE 
model with c = 0.8. In the upper panels, we show the cases of varying the neutrino mass. 
We choose — 0; OA, and 1.0 eV as examples; at the same time, other parameters are 

fixed (including Aeff = 3.046). We find that for the CMB power spectrum the only observed 
differences are for 2 < £ < 50 due to the late ISW effect from the neutrino background 
evolution. For the matter power spectrum, we find that massive neutrinos can suppress the 
amplitude of P{k) for a wide range of scales, but obviously the effect on small scales is more 
evident than on large scales. In the lower panels, we show the cases of varying the value 
of Aeff- We choose the examples of AAgfr = 0, 0.5, and 1.0; the neutrino mass is fixed to 
be ^ my = 0.06 eV. We find that for both CMB power spectrum and matter power 
spectrum P{k) the effect of dark radiation is more evident on small scales. 

Note that in this example we only show some unrealistic cases with one parameter 
varied and others fixed; in practice in order to fit to observational data, when one parameter 
is changed, other parameters would also change due to some compensation and degeneracy 
effects. Also, in fact, the CMB alone is not very powerful for constraining neutrino mass and 
dark radiation, and it should be used in combination with other measurements of expansion 
history and growth of structure. The matter power spectrum is useful for constraining 
neutrino mass but in practice it has limitations due to the bias and the fact that only 
intermediate region can be accurately measured. Thus, better choice is to use other LSS 
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data such as the measurements of weak leasing and redshift space distortions. 

In the next section, we will use the observational data to constrain the HDE model with 
massive neutrinos and/or dark radiation. We will use the Planck CMB data in combination 
with other geometric and structural growth measurements. For the geometric measurements, 
we use the BAO data, SN (JLA) data, and Hq measurement. For the measurements of growth 
of structure, we use the data of WL and RSD. 

3 Cosmological constraints 

We place constraints on the models of HDF with massive neutrinos and/or extra dark radia¬ 
tion by using the observational data. The conventions used in this paper are consistent with 
those adopted by the Planck team [7], i.e., those used in the camb Boltzmann code [48]. The 
base parameter set for the basic 7-parameter HDF model is: 

P = {ujb, ujc, lOO^MC, r, c. Us, ln(10^°As)}, 

where uJb = and Uc = are the baryon and cold dark matter densities today, 

respectively, 9mc is the approximation used in CosmoMC [49] to rs(z*)/d^(z*) (the angular size 
of the sound horizon at the time of last-scattering), r is the Thomson scattering optical depth 
due to reionization, c is the parameter governing the evolution of HDF, and Ug and Ag are 
the power-law spectral index and power amplitude of the primordial curvature perturbations, 
respectively. Flat priors for the base parameters are used. Note also that the prior ranges 
for the base parameters are chosen to be much wider than the posterior in order not to 
affect the results of parameter estimation. We use the CosmoMC package to infer the posterior 
probability distributions of parameters. 

3.1 Observational data 

Here, we describe the observational data sets used in this paper. ^ The data sets we use 
include the CMB, BAO, SN, i7o, WL, and RSD. 

The CMB data: We use the CMB TT angular power spectrum data from the 2013 
release of Planck [7], combined with the CMB large-scale TE and EE polarization power 
spectrum data form the 9-yr release of WMAP [6] . 

The BAO data: We use the BAO measurements from the 6dFGS {z = 0.1) [50], SDSS- 
DR7 {z = 0.35) [51], BOSS-DRll (z = 0.32 and 0.57) [52], and WiggleZ {z = 0.44, 0.60, and 
0.73) [53] surveys. This combination of BAO data has been used widely and proven to be in 
good agreement with the Planck CMB data. 

The SN data: We use the JLA compilation of the type la supernova observations con¬ 
taining 740 SN data [54]. 

The Hq measurement: We use the direct measurement of the Hubble constant from the 
HST observations, Hq = (73.8 ± 2.4) km s“^ Mpc“^ [55]. 

The WL data: We use the cosmic shear measurement of weak leasing from the CFHTLenS 
survey, cJ8(Dm/0.27)^'® = 0.79 ± 0.03 [56]. Moreover, since the CMB lensing reconstruction 

^In our description, there are lots of acronyms. For convenience, here we give a glossary of acronyms 
used in the description of observational data. WMAP: Wilkinson Microwave Anisotropy Probe; 6dFGS: Six- 
degree-Field Galaxy Survey; SDSS: Sloan Digital Sky Survey; DR: Data Release; BOSS: Baryon Oscillation 
Spectroscopic Survey; JLA: Joint Light-curve Analysis; HST: Hubble Space Telescope; CFHTLenS: Canada- 
France-Hawaii Telescope Lensing Survey; LRG: Luminous Red Galaxy; VIPERS: VIMOS Public Extragalactic 
Redshift Survey. 
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Figure 3. The CMB+BAO constraints on the models of HDE+^Wj^, HDE+A"eff, and 
HDE+y~] TTljy+Aeff ■ 


data directly probe the lensing power, thereby also sensitive to neutrino mass, we also use 
the CMB lensing power spectrum from the Planck mission [57]. 

The RSD data: We use the RSD measurements of f{z)as{z) from GdFGS {z = 0.067) 
[58], 2dFGRS {z = 0.17) [59], WiggleZ {z = 0.22, 0.41, 0.60, and 0.78) [60], SDSS LRG DR7 
{z = 0.25 and 0.37) [61], BOSS CMASS DRll {z = 0.57) [62], and VIPERS {z = 0.80) [63]. 

The CMB data contain both information of expansion and growth, so it is fairly im¬ 
portant for constraining the cosmological parameters. The BAO data probe the Hubble 
expansion rate and angular diameter distance at different redshifts, and have been proven 
to be in good agreement with the CMB data in a model-independent manner. Thus, we use 
the CMB-I-BAO as the basic data combination for constraining cosmological parameters in 
the HDE models. We consider the further geometric constraints by adding the SN (JLA) 
and Hq data in the cosmological fits. Eurthermore, we consider the constraints from growth 
of structure by adding the WL and RSD data in the fits. We consider three models in this 
paper, i.e., the HDE-|-^ mu model, the HDE-|-Aeff model, and the HDE-|-^ mu+Nes model. 
We will report the constraint results in the next subsection. 

3.2 Results and discussion 

We report the results for the HDE-|-^mi/ model (Model I), the HDE-|-AefT model (Model 
II), and the HDE-|-^ model (Model III) fitting to the three data combinations, 

i.e., the CMB-I-BAO (DI) combination, the CMB-|-BAO-|-SN-|-77o (D2) combination, and 
the CMB-|-BAO-|-SN-|-i7o+WL-|-RSD (D3) combination. Eor convenience, in the text of this 
subsection, we use Model I, Model II, and Model III to denote the three models, and use Dl, 
D2, and D3 to denote the three data combinations. 

The constraint results are shown in Eigs. 3-5 and Tables 1-3. Eig. 3 and Table 1 are 
responsible for the Dl constraints, Eig. 4 and Table 2 are for the D2 constraints, and Eig. 5 
and Table 3 are for the D3 constraints. In each figure, we show the two-dimensional posterior 
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Figure 4. The CMB+BAO+SN+i^o constraints on the models of HDE+^mjy, HDE+A^eff, and 
HDE+y^ mi/+A^eff ■ 
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Figure 5. The CMB+BAO+SN+i7o+WL+RSD constraints on the models of HDE+^m^, 
HDE+iVeff, and HDE+^ m^+A"eff ■ 


distribution contours (68% and 95% CL) in the Om-c, Qm-o's, Ngfi-Ho, and '^rriu- 

ag planes. Contours in green, orange, and blue stand for Model I, Model II, and Model 
III, respectively. One can easily compare different models under the same data combination 
in each individual figure, and it is also convenient to make a comparison for the constraint 
results of the same model with different data combinations by comparing the corresponding 
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0.982 ±0.014 

0.975 

0.981+°+^ 

c 

0.462 

0 c;+o.060 

U.010_Q 075 

0.571 

u.oyo_o.io3 

0.524 



0.033 

< 0.250 



0.007 

< 0.299 

A^eff 



3.45 

3.56 ±0.36 

3.54 

3.57 ±0.36 


0.245 

U.ZDi_0.016 

0.265 

n 071 +0.018 

U.Z ( 1_0.016 

0.255 

0.269 ±0.017 

as 

Ho 

0.897 

76.1 

0.857 ±0.036 
73.6^31 

0.863 

74.3 

f) «c;7+0-033 
n.OO/_0,039 

TA.Otl'r 

0.899 

76.4 

0 ft^O+0 037 

U.OOO_o.o40 

74.4tft 

-^max 


4906.71 


4906.56 


4906.47 


Table 1. Fitting results from the CMB+BAO data. We quote the ±1 ct errors, but for the neutrino 
mass '^rriu, we quote the 95% CL upper limits. Note that is in unit of eV and Hq is in unit 

of km s“^ Mpc“^. 




HDE±Em^ 

HDE±Afeff 

HDE±E ail, ± A^eff 

Parameter 

Best fit 

68% limits 

Best fit 

68% limits 

Best fit 

68% limits 


0.02228 

0.02234 ± 0.00025 

0.02292 

u.uzz 1 ^_o.00038 

0.02271 

0.02283 ± 0.0003 

1OO0MC 

T 

0.116 

1.04181 

0.105 

0.1151 ±0.0017 

i.UUiyO_g ggggg 

f) nqq+0 013 
u.uyy_g g 35 

0.1303 

1.04011 

0.114 

U.1ZD^_0.0072 

1 A^no+0.00067 

i.(J4U5_o. 00085 

U.lUO_o.oi9 

0.1256 

1.04059 

0.107 

0.1284;oo°g 

1.04062;°;" 

0 107+0-015 
U.iU/_Q 017 

Us 

0.9677 

0.9705 ± 0.0057 

1.008 

0.995 ±0.015 

0.992 

0.999 ±0.011 

c 

0.649 

n f;4Q+0-038 

U.D^.3_o.o44 

0.757 

n 71 r+0.053 

U. 1 -LO_o.o69 

0.713 

A 799+0.052 

U./ZZ_o 064 

Y.aiv 

0.0009 

< 0.133 



0.009 

< 0.207 

A^eff 



4.18 

3.841^2 

3.77 

q A 7 +O. 3 O 
'^•^‘-0.31 


0.2895 

0.2859 ± 0.0084 

0.2882 

0.289 ± 0.0083 

0.2871 

0.2891 ±0.0081 

as 

0.823 

0.808 ±0.021 

0.823 

n S91 +0.024 
U.OZi_Q Q26 

0.830 

A 099 + 0.026 

U.OZZ_o 023 

Ho 

69.1 

69.5l(i 

73.0 

72.1 ± 1.8 

71.9 

72.5 ± 1.4 

^71 /^niax 


5259.69 


5256.35 


5256.23 


Table 2. Fitting results from the CMB+BAO+SN+iLg data. We quote the ±lcr errors, but for the 
neutrino mass ^ we quote the 95% CL upper limits. Note that ^ rrii, is in unit of eV and iLg is 
in unit of km s“^ Mpc“^. 


panels in different figures. 

First, we discuss the constraint results from the D1 combination. From Fig. 3, one 
can see that there are strong parameter degeneracies in all the parameter planes. This is 
mainly due to the fact that using only CMB and BAO data cannot accurately measure these 
cosmological parameters in a dynamical dark energy model with neutrinos/dark radiation. 
The Planck data accurately measure seven acoustic peaks, and so the observed angular size of 
acoustic scale 0 * = TsIDa is determined to a very high precision (better than 0 . 1 % precision 
at la), implying tight constraints on some combinations of the cosmological parameters that 
determine and Da- The degeneracies in the Hm-c and Qrn~Ho planes can be viewed as 
mainly coming from the fact that the parameter combinations fitting to the Planck data must 
be constrained to be close to a surface of constant 0*. The BAO data have helped reduce 
the geometrical degeneracies to some extent, but not enough. Additional data are needed to 
further break the degeneracies, which will be seen in the following discussions. 
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HDE+2 


HDE±Neff 

HDE± JO nil, ± Aeff 

Parameter 

Best fit 

68% limits 

Best fit 

68% limits 

Best fit 

68% limits 


0.02225 

0.02249 ± 0.00024 

0.02276 

0.02292 ± 0.0003 

0.02297 

0.02292 ± 0.0003 


0.1145 

0.1137 ±0.0013 

0.122 

0.1234lHm5 

0.1236 

0.1239 ±0.0043 

lOO^MC 

T 

1.04165 

0.092 

1.04203 ± 0.00056 

U.lUZ_o 014 

1.04076 

0.115 

1.04098l°;ro^^ 

0 107+0014 

U.1UI_0.016 

1.04119 

0.109 

l.U^UyO_Q QQQgg 

0 111 +0-0^^ 

Us 

0.9700 

n Q7 Qc;+0-0056 

U.y i 00_Q QQ55 

0.992 

0.996 ±0.011 

0.996 

0.998 ±0.011 

c 

0.658 

0.670 ± 0.040 

0.728 

0 747+0-048 

u. < + <-0.060 

0.769 

n 7OQ+0.051 

U. ^ ZO_Q Qgg 


0.027 

< 0.186 



0.003 

< 0.343 

NeS 



3.68 

o 7r+0.28 
^ ^-0.32 

3.78 

3.79 ±0.30 


0.2860 

0.2869 ± 0.0080 

0.2846 

0.28791°;°°^^ 

0.2852 

n 0QQ'7+6-6U78 

QQg5 

CTS 

0.800 

0.792 ±0.014 

0.809 

0.802;nM 

0.810 

0.797t°;°[^ 

Ho 

69.21 

69.10t°|4 

71.5 

71.4 ± 1.3 

71.7 

7L7;]-i 

-^max 


5270.86 


5268.35 


5268.29 


Table 3. Fitting results from the CMB+BAO+SN+iJo+WL+RSD data. We quote the ilcr errors, 
but for the neutrino mass we quote the 95% CL upper limits. Note that is in unit of 

eV and Hq is in unit of km s“^ Mpc“^. 


We find that for all the three models, the degeneracy directions are the same. But 
we find that in the Llm-c and planes the contours of Model I (green) are somewhat 

lower, comparing to those of Model II (orange). Both neutrino mass and extra radiation 
could affect the shape of the CMB power spectra to some degree, but their effects are subtly 
different, as mentioned in the previous section. The effects on the background cosmology can 
be compensated by changes in, say, Hq (and c), to ensure the same observed acoustic peak 
scale 0*. Due to the subtle difference in the compensation effects, the contours of Model I 
and Model II exhibit some differences. Comparing the fittings of the two models, we find 
that Model II is better than Model I by = 0.3. Though Model III contains two extra 
parameters, ^ rrii, and N^s, it is only slightly better than Model II in the fit by = 0.18, 
implying that Agfr dominates over ^ niu in the fit and explaining why the contours of Model 
II and Model III are similar. 

In the case of D1 constraints, we find that c = O.SlS^^Qg^g for Model I, c = O.SOSlg'^gg 
for Model II, and c = 0.578 ^q5io for Model III. We thus find that lower values of c are derived 
once we only use the CMB and BAO data, consistent with the results of Ref. [40] . We can 
constrain Wgfr well but can only obtain an upper limit for ^ rriu. We find that ^ rriiy < 0.250 
eV (95% CL) for Model I and ^ < 0.299 eV (95% CL) for Model III; Nes = 3.56 ± 0.36 

for Model II and = 3.57 ± 0.36 for Model III. 

Second, we discuss the constraint results from the D2 combination. From Fig. 4, we 

find that all the contours are shrunk evidently, indicating that the constraint results are 

greatly tightened by adding more data of expansion history, especially the SN JLA data. It 
is clear to see that in this case the degeneracies are broken substantially due to the use of 
the SN data that have outstanding power for probing the properties of dark energy. In this 
fit, we find that Model II is much better than Model I by Ay^ = 6.68, but Model III is only 
slightly better than Model II by Ay^ = 0.24. This feature strongly indicates that in the fit 
Nes dominates over as demonstrated in the previous case. 

Under the constraints of D2 combination, for the HDE parameter, we find that c = 
0.643^0 Q 44 for Model I, c = 0.715lg ggg for Model II, and c = 0.722Y^'qq^ for Model III. So one 
can see that the values of c are enhanced once more geometric measurements data are added. 
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For the parameters of massive neutrinos and dark radiation, we find that '^rrii, < 0.133 eV 
(95% CL) for Model I and < 0-207 eV (95% CL) for Model III; A^eff = 3.84+°;^^ for 

Model II and A^eff = 3.97lg 3 ° for Model III. So, in this case, AA^eff > 0 is constrained to be 
at about the 3 (j level. 

Finally, we discuss the constraint results from the D3 combination. From Fig. 5, 
we find that the contours are further slightly shrunk by adding the measurements data of 
growth of structure, WL and RSD. Comparing Figs. 3-5, we notice that the inclusion of 
additional background data (SN+LIq) gives a visible improvement in the constraints, while 
the improvement is not evident when the growth data (WL+RSD) are further added. This 
is because we are considering a smooth dark energy model where the perturbations are 
suppressed on small scales. For the density perturbations in the HDE model, recall the 
examples shown in Fig. 1. Thus, the constraints from the structure growth probes are 
weaker. But we indeed find that in this case the fit values of erg ^-re further suppressed due 
to the consideration of the growth data. In this fit. Model II is better than Model I by 

Ax^ = 5.02, and Model III is better than Model II by Ay^ = 0.12. 

For this case, we find that c = 0.670 ± 0.040 for Model I, c = 0.747lQggQ for Model II, 
and c = 0.728^QQgg for Model III; we find that '^rriy < 0.186 eV (95% CL) for Model I and 
ruy < 0.343 eV (95% CL) for Model III; A"eff = 3 . 751 q 32 for Model II and Aeff = 3.79T0.30 
for Model III. In this case, AA^eff > 0 is still constrained to be at more than 2a level. Since the 

LSS data (WL and RSD) prefer lower ug, this leads to the upper limit on (comparing 

to the second case) being enhanced slightly. 

4 Conclusion 

In this paper, we investigated the models of HDE with massive neutrinos and/or dark ra¬ 
diation in detail. We calculated the background and perturbation evolutions in the HDE 
model. When the parameter c of HDE is less than 1, the EOS parameter w will evolve across 
the phantom divide w = — 1 , which usually leads to the perturbation divergence of dark 
energy. In order to overcome this instability difficulty, we employed the PPF approach in 
the perturbation calculations. We derived the evolutions of density perturbations of various 
components and metric fluctuations in the HDE model. We further discussed the impacts of 
massive neutrino and dark radiation on the CMB anisotropy power spectrum and the 
matter power spectrum P{k) in the HDE scenario. 

Furthermore, we constrained the models of HDE + HDE -|- Aeff; and HDE + 

Y + Aeff (denoted as Model I, Model H, and Model HI, respectively) by using the data 
combinations of CMB-^BAO, CMB+BAO+SN-bi7o, and CMB-FBAO-KSN+iLo+WL-^RSD 
(denoted as Dl, D2, and D3 combinations, respectively). With the latest measurements of 
expansion history and growth of structure, we found that the HDE models can be tightly 
constrained. 

When we considered the Dl combination, we found that there are strong degeneracies 
among the cosmological parameters. Since the Planck data accurately determine the acoustic 
scale O-t to a very high precision, some parameter combinations that determine and Da 
are tightly constrained, but this also leads to some degeneracies due to the fact that the 
parameter combinations must be constrained to be close to a surface of constant 0*. Though 
with the help of BAO data, some degeneracies have been reduced to some extent, this is 
not enough, and more additional data are needed to further break the degeneracies. Thus 
we further considered the additional background expansion data, SN (JLA) and Hq. Under 
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the constraints of D2, we found that the degeneracies are broken well. This indicates that 
for constraining a dynamical dark energy model (with neutrinos/dark radiation), the back¬ 
ground data, especially the SN data, play a crucial role. Furthermore, we also considered the 
structure growth data, WL and RSD. We found that in the D3 fit the constraints are only 
slightly improved comparing to the D2 case, showing that comparing with the expansion 
probes, the constraints from growth probes are weaker. This is because we are considering 
a smooth dark energy model where the perturbations are suppressed on small scales. In 
addition, in all the cases, we found that in our cosmological fits Model II is much better than 
Model I, while Model III is only slightly better than Model II (though the former with one 
more parameter). This implies that dominates over ™ 

Using the CMB-|-BAO-|-SN-|-FIo+WL-|-RSD data, we obtained the results: c = 0.670 ± 
0.040 and < 0.186 eV (95% CL) for the HDE-|-^m,y model, c = 0.747lgQQQ and 

Aeff = 3.75^0 32 for the HDE-|-Aefj model, and c = 0.728^0000, ^ rriy < 0.343 eV (95% CL), 
and A^eff = 3.79 ± 0.30 for the HDE-|-^ mv+N(,Q model. 
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